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Key Points 


è The global response of temperature of cold-point-mesopause to solar activity is 4.89+0.67 
K/100SFU 

+ Temperature of cold-point-mesopause and its solar response increases gradually with 
increasing latitude 

+ Latitudinal distribution of correlation coefficients of annual-mean series between solar flux 


and temperature of mesopause takes on M-shape 


Abstract 

The thermal structure and energy balance of upper atmosphere are dominated by solar 
activity. The response of Cold-Point-Mesopause(CPM) to solar activity is an important form. 
This article presents the response of Temperature-of-CPM(T-CPM) to solar activity using 
fourteen-year SABER data series over 80°S—80°N regions. These regions are divided into 16 
latitude zones with 10° interval, and the spatial areas of 80°S—80°N, 180°W—180°E are divided 
into 96 lattices with 10°(latitude)x60°(longitude) grid. The annual-mean values of T-CPM and 
Fio7 are calculated. The least squares regression method and correlation analysis are applied to 
these annual-mean series. First, the results show that the global T-CPM is significantly correlated 
to solar activity at the 0.05 level of significance with correlation coefficient of 0.90. The global 
solar response of T-CPM is 4.89+0.67K/100 Solar-Flux-Units(SFU). Then, for each latitude 
zone, the solar response of T-CPM and its fluctuation are obtained. The solar response of T-CPM 
becomes stronger with increasing latitude. The fluctuation ranges of solar response at middle 
latitude regions are smaller than those of equator and high latitude regions, and the global 
distribution takes on W-shape. The co-relationship analysis shows that the T-CPM is 
significantly correlated to solar activity at the 0.05 level of significance for each latitude zone. 
The correlation coefficients at middle latitude regions are higher than those of equator and high 
latitude regions, and the global distribution takes on M-shape. At last, for each grid cell, the 


response of T-CPM to solar activity and their correlation coefficient are presented. 
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1 Introduction 

The energy balance of upper atmosphere of earth is greatly affected by solar radiation. The 
most energetic solar photons at ultraviolet (UV) and shorted wavelengths are removed at various 
levels in the atmosphere, and the earth’s surface re-emits the absorbed energy at long 
wavelength. So many energy source and sinks processes are occurred at upper atmosphere, and 
its thermal structure is greatly impacted by the abundance of solar energy striking the Earth from 
a constantly changing sun. Variations arising on decadal and even longer timescales may play a 
significant role in long-term trend estimates. One of the major sources of decadal variability in 
the upper atmosphere is the 11-year solar activity cycle. Electromagnetic radiation from the Sun 
is not constant and varies mainly at shorter ultraviolet wavelengths on different timescales. 
Incoming solar radiation provides the external forcing for the Earth-atmosphere system. While 
the total solar flux is quite constant, the UV spectral irradiance on the timescale of the 27-d and 
11-year solar cycles exhibits the largest changes, up to a factor of 2 over a solar cycle for the 
solar Lyman alpha flux [Beig et al., 2008]. The thermal structure and energy balance of upper 
atmosphere are dominated by solar activity cycle. The response of the cold-point-mesopause to 
solar activity variations is an important form. 

The Cold-Point-Mesopause (CPM), defined as the region where kinetic temperature is 
minimum, separates the mesosphere from the thermosphere, in which the transport and exchange 
of energy occur through subtle and complex processes [Beig et al., 2003]. Temperature is the 
most important parameter for characterizing this transition from the mesosphere to the 
thermosphere, and the variation of temperature may be related with the solar activities. 
Temperature-of-CPM (T-CPM) and Height-of-CPM (H-CPM) are established by radiative, 
chemical and dynamical processes, and also display large variability due to internal gravity 
waves and tides [Xu et al., 2007; John et al., 2011]. 

Solar energy is the original source of all these complex atmospheric processes. Hence, the 
investigations of the changes in solar UV spectral irradiance on timescales of the 27-day and 11- 
year solar cycles have been studied by many researchers. She et al. [2009] investigated solar 
cycle effect and temperature trends in mesopause region on the basis of 894 nights of Na lidar 
observation obtained during the period 1990-2007 over Fort Collins (41°N, 105°W). Offermann 
et al. [2010] analyzed the response of temperature to solar activity about ~87 km on the basis of 
OH airglow measurements obtained during the period 1988—2008 over Wuppertal station (51°N, 
7°E). Ammosova et al. [2012] investigated temperature data on the basis of spectral observations 


of the mesopause emissions on Maymaga station (63°N, 129°E) during the 2000-2011 years 
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seasonal dependence of the response received radiation temperatures of OH (6,2) (height of the 
radiation 87 km) from the solar activity. Savigny et al. [2012] studied the sensitivity of equatorial 
mesopause temperatures to the 27-day solar forcing using night-time observations of OH(3-1) 
rotational temperatures with SCIAMACHY (SCanning Imaging Absorption spectroMeter for 
Atmospheric CHartographY) on Envisat. Ammosov et al. [2014] analyzed the measurements of 
rotational temperature of the hydroxyl near the mesopause region at Maimaga station (63.04°N, 
129.51°E) for the period from August 1999 to March 2013 in search for solar activity effects and 
temperature trends. However, the above literatures mainly reported some separate local 
mesopause region through analyzing temperature data obtained by ground-based remote sensing. 
The ground-based remote sensing, such as lidar measurements or OH airglow measurements, can 
only observe a small portion of the upper atmosphere located over an observation site. So, the 
ground-based remote sensing was based on separate local sites with a limited number of 
observation days, and failed to give continuous observations on global scale. 

Satellite remote sensing, such as HALOE observations or SABER observations, is the most 
efficient method to obtain the atmosphere parameters for investigating the global scale effects of 
solar variability on temperature in the upper atmosphere. Remsberg [2009] reported solar cycle 
effects in temperature versus altitude from the Halogen Occultation Experiment (HALOE) for 
the mesosphere. Oindrila and Sridharan [2014] investigated the long-term variabilities and trends 
of middle atmospheric temperature and their responses towards quasi-biennial oscillation(QBO), 
solar cycle and EINifio-southern oscillation(ENSO), extending 20—100 km altitude and 10—15°N 
latitude using SABER observations for the years 2002-2012. Forbes et al. [2014] investigated a 
response of the middle atmosphere at various heights on change of solar activity using SABER 
temperature data during 2002-2013 and covering 20-110 km altitude. Ramesh et al. [2015] 
studied the solar cycle variabilities of tropical (10°N—15°N) mesopause temperature structure 
through analyzing SABER observations during 2002-2012. Detailed reviews of this subject were 
made by Beig et al. [2008] and Beig [2011]. However, the global distribution of the response of 
T-CPM (at mesopause actually) to solar activity, which was obtained using same data source 
with coverage of at least one complete solar cycle, was not reported. In our work, the close 
connection of T-CPM with the level of solar activity is statistically presented on the basis of 
large volume of data, and the nonuniform latitudinal and longitudinal distribution of the response 
of T-CPM to solar activity is shown. 

The Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) is one 
of four instruments on NASA's TIMED (Thermosphere Ionosphere Mesosphere Energetics 
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Dynamics) satellite launched in December 2001. The primary goal of the SABER experiment is 
to provide the data needed to advance our understanding of the fundamental processes governing 
the energetics, chemistry, dynamics, and transport in the mesosphere and lower thermosphere. 
The latitudinal coverage of SABER observations is either 83°N to 52°S or 83°S to 52°N 
depending on the yaw cycles, and it makes ~15 orbits per day with a period of ~1.6 hour. The 
SABER scans the atmosphere from the troposphere of the Earth up to lower thermosphere once 
in every ~58s, and obtains vertical profiles of kinetic temperature, geopotential height and other 
parameters [Russell et al., 1999; Mertens et al., 2001,2003; Remsberg et al, 2003]. The SABER 
observations offer sufficiently good quality dataset, and have larger scope and much longer 
records than OH airglow measurements and lidar measurements. Thus, SABER observations are 
optimal data source for investigating the global scale long-term temperature response to solar 
activity. According to the characteristics of SABER observations, for the convenience of division 
of latitude zones and analysis, we will present the response of T-CPM to solar activity and their 
correlation coefficient globally from 80°S to 80°N based on 2002-2015 SABER temperature 
products in this paper. The solar activities are presented with the 10.7cm solar radio flux (Fio7) 
dataset. The results are restricted to nighttime situation to minimize the impact of fluctuations 
caused by the propagation of internal gravity waves and tides. In section 2, the data and 
methodologies used in the investigation will be introduced. The results of the response of T- 
CPM to solar activity will be shown in section 3, and a brief discussion and conclusions will be 


given in last section. 


2 Data and Methodologies 
2.1 Data source 

The atmospheric temperature profiles we used in this investigation are SABER version 
2.0(V2.0) Level 2A dataset, which lasts from 2002.1.25 up to now. The V2.0 products were 
available to public in November 2012. Compared to prior versions, the V2.0 production retrieval 
algorithm has many substantial improvements. These improvements significantly impact all 
parameters. More details on retrieval algorithm and data quality assessment can be obtained from 


Mertens et al.[2003], Remsberg et al. [2003, 2008] and the website: http://saber.gats- 


inc.com/V107vsV2.php. 


The 10.7 cm solar radio flux (F107), which varies with solar cycle, is often taken as an index 
of solar activity [e.g. Gray et al., 2010; Beig, 2011; Ammosova et al., 2012; Ammosov et al., 
2014; Ramesh et al., 2015]. The daily data of Fio7 adjusted for an arbitrary unit (1 A.U.) is used 


to analyze the response of T-CPM to solar activity. 


2.2 The methodology for CPM data quality control processing 

Through processing 14 years (14yr) temperature profiles data of SABER observations from 
2002 to 2015, the T-CPM is firstly identified for each temperature profile based on the definition 
of CPM. The H-CPM is corresponding height for the T-CPM. Owing to some incorrect CPMs in 
the temperature profiles, it is necessary to preprocess the profiles by carrying out quality control 
for T-CPM. For all the 4,495,286 temperature profiles with any valid nighttime event reports in 
level 2A dataset of the SABER version 2.0 during 2002 to 2015, the following three steps are 
carried out for each temperature profile. Firstly, the temperature profile with the maximum 
height lower than 90 km, which cannot give an effective judge for T-CPM and H-CPM, is 
ignored. Secondly, the temperature profile is ignored when height of the temperature profile does 
not rise above the H-CPM for it is an ineffective record, as the circle symbol shown in Fig. la. 
Finally, the circle symbol shown in Fig.lb is considered as an effective CPM, and the 
corresponding effective temperature, height, date, latitude and longitude are saved for the 


following analysis. 
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Figure 1. Examples of the CPM data quality control process. (a) The temperature profile is the 100th event in 
75361 orbit at 2015-11-4 night; (b) The temperature profile is the 40th event in 76214 orbit at 2015-12-31 night 


Fig.2a shows the distribution of amount of global available temperature profiles for each 
year from 2002 to 2015 in level 2A dataset of the SABER version 2.0. The original data are 
shown in hollow hexagon, and the effective data after CPM quality control are shown in solid 
circle. After the quality control of CPM, the total amount of global available temperature profiles 
of T-CPM decreases from the original 4,495,286 to 4,492,664. A total of 2622 temperature 
profiles are removed in the outside of CPM quality control for T-CPM, and approximately 99.94 
percent of 4,495,286 profiles are effective values. Also found in Fig.2a, due to SABER 
observations lack of data during 2002.1.1 to 2002.1.24, the amount of temperature profiles of T- 
CPM in 2002 is about 1/14 less than that of other years. In order to match with the satellite data 
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series, the daily data of Fio7 during 2002.1.1 to 2002.1.24 are also removed out. 

Fig.2b shows the distribution of 4,492,664 temperature profiles of T-CPM for each integer 
latitude interval between the latitudes of about 83°S to 83°N. From the solid circle symbol shown 
in Fig.2b, we can find that, between 52°S and 52°N, the amount of profiles of T-CPM are more 
than 25000 for each integer latitude interval, and the amount of profiles of T-CPM are more than 
10000 between 52°N and 83°N or between 52°S and 83°S. The uneven distribution of 
temperature profiles accords with the latitudinal coverage of SABER observations which is 
either 83°N to 52°S or 83°S to 52°N depending on the yaw cycles. Latitudinal coverage varies 
with time because SABER performs a yaw maneuver approximately every 60 days, so only the 
latitude range from 52°S to 52°N is continuously covered by SABER observations. Based on 
above analysis, in order to avoid the yaw cycles impact on the solar response of T-CPM, we only 
analyze the annual-mean response of T-CPM to solar activity and their correlation coefficient, 
and the global annual-mean of T-CPM is calculated by averaging latitude zones mean values in 


this study. 
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Figure 2. Distribution of global available profiles of T-CPM. (a) Amount of total temperature profiles of T- 
CPM for each year before CPM quality control (hollow hexagon) and after CPM quality control (solid circle); (b) 
Distribution of 4,492,664 profiles of T-CPM for each integer latitude interval after CPM quality control 


At last, during the period from 2002 to 2015, the total number of 4,302,764 temperature 
profiles are adopted in the following investigations over 80°S—80°N latitude region, which is 
about 45 times of the amount of HALOE temperature profiles reported in the literature 
[Remsberg, E. E., 2009]. 

2.3 The methodologies for CPM data analyzing 

We divide the latitude regions from 80°S to 80°N into 16 latitude zones with the interval of 
10°. The annual-mean of T-CPM for each latitude zone is calculated by the T-CPM of all 
temperature profiles within the latitude zone. For the investigations of the spatial (latitude- 
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longitude) distribution of response of T-CPM to solar activity and their correlation coefficient, 
we divide the areas of 80°S—80°N, 180°W—180°E into 96 lattices using 
10°(latitude)x60°(longitude) grid to ensure each grid cell contains a sufficient number of CPM 
samples. The annual-mean of T-CPM for each grid cell is calculated by the T-CPM of all 
temperature profiles within the cell. 

In this study, the global annual-mean of T-CPM is calculated by the corresponding average 
of 16 latitude zones, and the annual-mean of Fio7 is calculated by the daily data of Fio7. Using 
least squares method and linear regression analysis, we obtain the global solar response of T- 
CPM and the fluctuation range of solar response. Through correlation analysis, we obtain the 
global correlation coefficient of annual-mean series values between Fi); and T-CPM, and the 
global significance of correlation coefficient is tested with 12 degrees of freedom T-test. Then, 
for each latitude zone and each grid cell, the solar response of T-CPM and its fluctuation are 
analyzed by the least squares regression method, and the corresponding significance of 
correlation coefficient is also tested with 12 degrees of freedom T-test. The spatial distribution 


characteristics of solar response of T-CPM and their solar cycle dependence are obtained. 


3. Results 


3.1. Global response of T-CPM to solar activity and their correlation coefficient 

In order to reduce the noises and increase the statistical significance, the global annual- 
mean of T-CPM is calculated by the corresponding average of 16 latitude zones. By using least 
squares method and linear regression analysis, the scatter plot of T-CPM and Fio7 is shown in 
Fig.3a, the global solar response of T-CPM and the fluctuation range of solar response are 
obtained. From Fig.3a, we can find that the T-CPM varies almost linearly with the solar flux 
Fio7. The global solar response of T-CPM is 4.89+0.67 K/100SFU. 

Through the calculation of correlation coefficient, we analyze the dependence relation of 
fourteen-year time series average between T-CPM and Fio7, as shown in Fig.3b, and the 
correlation coefficient is 0.90. Since the probability density function of sample correlation 
coefficient obeys the T-distribution [Mises et al., 1964], the significance of correlation 
coefficient will be tested with 12 degrees of freedom T-test in this paper. By calculation, the 
critical value of correlation coefficient is 0.53 at the 0.05 level of significance. The co- 
relationship analysis shows that the global annual-mean of T-CPM is significantly correlated to 
the solar activity at the 0.05 level of significance during 2002-2015. Solar variations are the 


dominant influence on the upper atmospheric temperature structure. 
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Figure 3. The scatter plot and the annual-mean of T-CPM and solar flux index during 2002-2015. (a) The 
scatter plot of T-CPM and Fio7; (b) The annual-mean of T-CPM and Fio7 


Fig.3b clearly shows that there is a strong dependence of the T-CPM on the changes of solar 
activity. The T-CPM decreases from 2002 to the bottom at 2008, and then gradually increases to 
the top. The changes of T-CPM are closely with the changes of solar activity. It is commonly 
known that the solar activity has 11-year solar activity cycles. The annual-mean of T-CPM 
coincides with the solar 1 1-year-cycle changes. 

One important mechanism of this dependency is that CO: cooling is dominant during solar 
minimum, while NO cooling becomes more important during maximum [Qian et al. 2006]. The 
solar cycle dependence is due to the relative importance of CO, infrared cooling and NO infrared 
cooling. NO density at solar maximum is about three times of that at solar minimum. 
Consequently, CO, cooling is relatively less important at solar maximum [Qian et al. 2011]. 
SABER measurements revealed a decrease of the global NO cooling rate from 2002 to 2009 by 
nearly an order of magnitude, while the global CO, cooling rate decreased only by 35% 
[Mlynczak et al., 2010]. 

The T-CPM is determined by many factors under the radiative, chemical and dynamical 
processes. Because the sun is the original source of all these complex atmospheric processes, it is 
not surprising that the changes of the solar energy result in the variation of these processes, and 
make the changes of T-CPM. We find the changes of T-CPM depend closely on the changes of 
solar activities from the SABER observations over the global region. 

With the increasing of solar flux index Fio7, the solar UV photons increase, the production 
of ion pairs by photoionization(N2,O2) and the production of atoms by molecular(O2) 
photodissociation increase, thus, the abundance of O atoms vary with solar activities and in turn 
induce the changes in the thermal structure of the middle atmosphere thereby affecting the 


temperature-dependent reaction rates, long- and short-wave radiative heating. In the lower 
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thermosphere, thermal collisions with O atoms provide the principal excitation process, so the 
temperature varies with the variation of solar activity. These may explain the reasons that the T- 
CPM varies with the solar activities. 

The global annual-mean of H-CPM is also calculated by the corresponding average of 16 
latitude zones, as shown in Fig.4. The annual-mean of H-CPM serves as the location reference of 


solar response of T-CPM. The variation of annual-mean of H-CPM is small, it is around 


95.4+0.4km, it has negative dependence with the solar cycle with the correlation coefficient of 
-0.78. 
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Figure 4. The annual-mean of H-CPM during 2002-2015 


3.2. Latitudinal distribution of the response of T-CPM to solar activity and their 
correlation coefficient 


In order to study the latitudinal distribution of the response of T-CPM to solar activity and 
their correlation coefficient, the annual-mean of T-CPM for each latitude zone is calculated by 
the T-CPM of all temperature profiles within the latitude zone. By using least squares method 
and linear regression analysis, the latitudinal distribution of the solar response of T-CPM and the 
fluctuation range of solar response are obtained. Fig.5a shows the solar response of T-CPM and 
the fluctuation range of solar response with every 10° interval. From Fig.5a, we can see that the 


latitudinal distribution characteristics of the solar response of T-CPM are remarkable. 
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Figure 5. Latitudinal distribution of solar response and solar cycle dependence with every 10° interval. (a) 
Latitudinal distribution of solar response of T-CPM, and the vertical bars show the fluctuation range of solar 
response; (b) Latitudinal distribution of correlation coefficients between T-CPM and F107 


In the northern hemisphere, the solar response of T-CPM increases gradually from low 
latitude zone to high latitude zone with the solar response values vary from the minimum 
2.78+40.69 K/100SFU in 0°N—10°N latitude zone to the maximum 7.21+1.85 K/100SFU in 
70°N—80°N latitude zone. The fluctuation range of solar response increases gradually from 
20°N—30°N latitude zone to high latitude zone with the fluctuation values vary from +0.35 
K/1OOSFU in 20°N—30°N latitude zone to +1.85 K/100SFU in 70°N—80°N latitude zone, and 
in equator region, the fluctuation range of solar response decreases gradually from 0°N—10°N 
latitude zone with +0.69 K/100SFU to 20°N—30°N latitude zone with +0.35 K/100SFU. The 
smallest fluctuation of solar response with +0.35 K/100SFU occurs in 20°N—30°N latitude 
zone, and the largest fluctuation of solar response with +1.85 K/100SFU occurs in 70°N—80°N 
latitude zone. 

In the southern hemisphere, except for 70°S—80°S latitude zone with 6.13+1.21 
K/100SFU, the solar response of T-CPM increases gradually from low latitude zone to high 
latitude zone with the solar response values vary from 2.82+0.73 K/100SFU in 0°S—10°S 
latitude zone to 6.35+1.16 K/100SFU in 60°S—70°S latitude zone. The minimum solar response 
with 2.82+0.73 K/100SFU occurs in equator latitude region (0°S—10°S), and the maximum 
solar response with 6.35+1.16 K/100SFU occurs in 60°S—70°S latitude zone. Except for 0°S— 
10°S latitude zone with +0.73 K/100SFU, the fluctuation range of solar response increases 
gradually from low latitude zone to high latitude zone with the fluctuation values vary from 
+0.42 K/100SFU in 10°S—20°S latitude zone to +1.21 K/100SFU in 70°S—80°S latitude zone. 
The smallest fluctuation of solar response with +0.42 K/100SFU occurs in 10°S—20°S latitude 
zone, and the largest fluctuation of solar response with +1.21 K/100SFU occurs in 70°S—80°S 


latitude zone. Based on above analysis and combined with Fig.5a, it can be concluded that the 
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fluctuation ranges of the solar response at middle latitude regions are remarkable smaller than 
those of equator and high latitude regions, and the solar response fluctuation with latitude at the 
global scale takes on W-shape. The phenomenon was not reported before. 

By using correlation analysis, we obtain the correlation coefficients of annual-mean series 
between Fio7 and T-CPM for 16 latitude zones, as shown in Fig.5b. The correlation coefficients 
also show obviously latitudinal distribution characteristics. The correlation coefficients at middle 
latitude regions are higher than those of equator and high latitude regions, and the correlation 
coefficients with latitude at the global scale take on M-shape. The phenomenon was also not 
reported before. Due to the correlation coefficient of each latitude zone greater than 0.53, the co- 
relationship analysis shows that the annual-mean of T-CPM is significantly correlated to the 
solar activity at the 0.05 level of significance for each latitude zone. 

Based on above analysis and combined with Fig.5, it can be concluded that the correlation 
coefficient depends on the corresponding fluctuation range of solar response. The lower the 
correlation coefficient, the larger the fluctuation amplitude. The phenomenon that the larger 
fluctuation and lower correlation coefficient at the equator and high latitude zones might display 
that the variation of temperature at the mesopause atmosphere is more active at the equator and 
high latitude zones. 

To interpret the reliability of above results, we count the number of CPM samples at each 
latitude zone. Fig.6 shows the distribution of global available CPM samples for each latitude 
zone. From Fig.6, we can see that each latitude zone contains a sufficient number of CPM 
samples, and the amount of CPM samples is more than 10° for each latitude zone. Thus, the 
number of CPM samples is enough large to account the statistics credible for the solar response 


of T-CPM over the global zones. 
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Figure 6. Latitudinal distribution of 4,302,764 temperature profiles with every 10° interval 
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Fig.7a shows the latitudinal distribution of 14yr average of T-CPM for each latitude zone, 
and the 14yr average of T-CPM is calculated by the T-CPM of all temperature profiles within the 
latitude zone during 2002-2015, and the corresponding H-CPM is shown in Fig.7b. The 14yr 
average of H-CPM serves as the location reference of solar response of T-CPM for each latitude 
zone. From Fig.7a, we can see that the latitudinal distribution characteristics of T-CPM are 
obviously, and the T-CPM increases gradually from the equator to both sides to high latitude 
zones. It is known that the behavior of a seasonal mesopause of high latitudes has a 2-level 
structure: mesopause drops to ~85 km and cools down to 130-140 K in summer. Due to the 
latitude zones contain uneven amount of CPM samples shown in Fig.6, and there is no CPM 
samples for some nights at polar latitude zones in summer, the 14yr average of T-CPM and H- 
CPM may be lower for high latitude zones above 60 degrees. 
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Figure 7. Latitudinal distribution of 14yr-average for T-CPM and H-CPM with every 10° interval. (a) Latitudinal 
distribution of 14yr average of T-CPM; (b) Latitudinal distribution of 14yr average of H-CPM 


3.3. Spatial distribution of the response of T-CPM to solar activity and their correlation 


coefficient 

In this section, we will focus on the spatial (latitude-longitude) distribution of the response 
of T-CPM to solar activity and their correlation coefficient for each grid cell. The annual-mean 
of T-CPM for each grid cell is calculated by the T-CPM of all temperature profiles within the 
grid cell. Using least squares regression method, the solar response of T-CPM is obtained for 
each grid cell. Through correlation analysis, we obtain the correlation coefficient of annual-mean 
series between F107 and T-CPM, and then the significance of correlation coefficient is tested with 
12 degrees of freedom T-test for each grid cell. Fig.8a shows the spatial distribution of solar 
response of T-CPM, and Fig.8b shows the spatial distribution of correlation coefficients of 


annual-mean series between T-CPM and F107 during 2002-2015. 


13 


(a) The Solar Response of T-CPM Unit: K/100SFU (b) Correlation Coefficient 
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Figure 8. Maps of the solar response and solar cycle dependence. (a) Spatial distribution of the solar response of T- 
CPM, and contour interval is 1 K/100SFU; (b) Spatial distribution of the correlation coefficients between T-CPM 
and Fjo.7, and contour interval is 0.05 above 0.53 


As shown in Fig.8a, the solar response of T-CPM shows obviously spatial distribution 
characteristics, and the variation amplitude of the solar response of T-CPM vary from 2.12 
K/100SFU to 8.18 K/100SFU. The minimum solar response of T-CPM with 2.12 K/100SFU 
occurs in the grid cell which is located at 0°N—10°N,0°E—60°E region, as the black star 
symbol shown. The maximum solar response of T-CPM with 8.18 K/100SFU occurs in the grid 
cell which is located at 70°N—80°N,60°W—120°W region, as the black triangle symbol shown. 
Under the same longitude zone with every 60° interval, on the whole, the change gradient of the 
solar response of T-CPM increases gradually from the equator to both sides to high latitude 
regions except for some individual grid cells. 

As shown in Fig.8b, the correlation coefficients also show obviously spatial distribution 
characteristics, and the variation amplitude of correlation coefficients between T-CPM and Fio7 
vary from 0.51 to 0.95. The minimum solar cycle dependence with the correlation coefficient of 
0.51 occurs in the grid cell which is located at 0°N—10°N,0°E—60°E region, as the black star 
symbol shown. The maximum solar cycle dependence with the correlation coefficient of 0.95 
occurs in the grid cell which is located at 30°N—40°N, 0°E—60°E region, as the black triangle 
symbol shown. For each grid, the significance of correlation coefficient is tested with 12 degrees 
of freedom T-test. Only one grid cell region, as the blue color shown in map, is not significantly 
correlated to the solar activity at the 0.05 level of significance. The remaining ninety-five grid 
cells are significantly correlated to the solar activity at the 0.05 level of significance. 

Also, both the solar response of T-CPM and their correlation coefficient show a little 
longitudinal dependence, as can be seen in Fig.8. 

Fig.9a shows the spatial distribution of 14yr average of T-CPM for each grid cell, and the 
14yr average of T-CPM is calculated by the T-CPM of all temperature profiles within the grid 
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cell during 2002-2015, and the corresponding H-CPM is shown in Fig.9b. The 14yr average of 
H-CPM serves as the location reference of solar response of T-CPM for each grid cell. Same as 
Fig.7a, the T-CPM shows the remarkable latitudinal dependence, and has little longitudinal 


dependence. 
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Figure 9. Maps of 14yr-average for T-CPM and H-CPM. (a) Spatial distribution of 14yr average of T-CPM, and 
contour interval is 2K; (b) Spatial distribution of 14yr average of H-CPM, and contour interval is 0.5km 


4. Discussion and Conclusions 


4.1. Discussion 

A wide scatter of literatures have investigated the solar activity influence on upper 
atmosphere temperature [Remsberg, 2009; She et al., 2009; Offermann et al., 2010; Ammosova 
et al., 2012; Savigny et al., 2012; Ammosov et al., 2014; Forbes et al., 2014; Oindrila and 
Sridharan, 2014; Ramesh et al., 2015]. Beig et al. [2008] and Beig [2011] provided a compilation 
of the existing studies quantifying the sensitivity of middle atmospheric temperatures to solar 
activity in terms of the 11-year solar activity cycle, the majority estimations lie in the range from 
0 to 10 K/100 SFU in mesopause region. We process a large array of data, and the results show 
that there is a strong dependence of the temperature of cold-point-mesopause and its latitudinal 
distribution from changes in solar activity. Our results, for northern and southern hemisphere that 
refer to heights around 95 km, coincide with most of the conclusions obtained by different 
observations separated by several decades. However, the magnitude of solar response reported 
by different authors is found to be a little different compared to our results. Some of the 
differences in temperature response to solar activity in the mesopause region are understandable. 
The differences may be attributed to the fact that ground-based remote sensing (such as lidar 
measurements and OH airglow measurements) results were averages of measurements obtained 
over several nighttime hours, whereas the HALOE results were strictly based on the data 


obtained for sunrise and sunset time only, and the HALOE observations were intermittent, as 
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shown from the website: http://haloe.gats-inc.com/coverage/index.php. The differences, the 
magnitude of solar response of T-CPM reported by other researchers compared to our results, 
may be due to the height difference since the solar temperature response increases with altitude 
from the middle mesosphere into the lower thermosphere [Forbes et al., 2014]. The differences 
may also be due to decadal or solar cycle effect in the tidal forcing at tropics, where tidal 
amplitudes are large [Beig, 2011]. 

Due to no more results in previous literatures, which were obtained using same data source 
with long-term data series, were available for the global scale mesopause region, the major 
challenge is in the interpretation of the latitudinal or longitudinal variability of the response of T- 
CPM to solar activity. Our results show that the spatial distribution characteristics of the solar 
response are obviously. The spatial changes of the response of T-CPM to solar activity may be 
mainly caused by the changes in the spatial distribution of chemically active gases (such as 
ozone, NO, CO», and O) and by changes in Solar UV irradiation. The spatial changes may also 
be related with the solar radiation intension or solar incidence angle in the difference mesopause 
region, as for the other reasons for the spatial changes are worthy of further investigations. 

4.2. Conclusions 

The first global distribution of the response of T-CPM to solar activity and their correlation 
coefficient are presented through analyzing temperature profiles of the level 2A dataset in the 
SABER version 2.0 from 2002 to 2015. The analyses include the latitudinal distribution, spatial 
distribution, the solar response of T-CPM and its fluctuation, and solar cycle dependence and its 
significance of correlation coefficient. The results show that there is a close connection of T- 
CPM (at mesopause actually) with the level of solar activity, and there exists a nonuniform 
latitudinal and longitudinal distribution of the response of T-CPM to solar activity on the basis of 
large volume of data. 

The 4,302,764 temperature profiles are applied for this period. By using the least squares 
method and linear regression analysis, the solar response of T-CPM and its fluctuation range are 
presented. The global solar response of T-CPM is 4.89+0.67 K/100SFU. The global T-CPM 
clearly follows the solar cycle, and the corresponding significance correlation coefficient is 0.90. 

Through the analysis for the latitudinal distribution of the response of T-CPM to solar 
activity and its fluctuation, it is found that the latitudinal distribution characteristics of the solar 
response are obviously. The solar response of T-CPM increases gradually from the equator to 
both sides to high latitude zones except for one latitude zone(70°S—80°S), and the fluctuation 


range of solar response increases gradually from low latitude zones to high latitude zones except 
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for three latitude zones in tropic region over 10°S—20°N. The co-relationship analysis shows 
that the annual-mean of T-CPM is significantly correlated to the solar activity at the 0.05 level of 
significance for each latitude zone. 

The following phenomena were not reported in previous literatures: the fluctuation ranges 
of solar response at middle latitude regions are remarkable smaller than those of equator and 
high latitude regions, and the solar response fluctuation with latitude at the global scale takes on 
W-shape. The correlation coefficients at middle latitude regions are higher than those of equator 
and high latitude regions, and the correlation coefficients with latitude at the global scale take on 
M-shape. The correlation coefficient and the corresponding fluctuation range of solar response 
are related values, and the lower the correlation coefficient, the larger the amplitude. 

Through the analysis for the spatial distribution of the response of T-CPM to solar activity 
and their correlation coefficient for each grid cell, the nonuniform latitudinal and longitudinal 
distribution of the solar response of T-CPM is shown. The values of solar response of T-CPM 
vary from 2.12 K/100SFU at the equator region to 8.18 K/100SFU at high latitude region. The 
minimum solar response of T-CPM with 2.12 K/100SFU occurs in the grid cell which is located 
at 0°N—10°N,0°E—60°E region. The maximum solar response of T-CPM with 8.18 K/100SFU 
occurs in the grid cell which is located at 70°N—80°N, 60°W—120°W region. For each grid 
cell, the significance of correlation coefficient is tested with 12 degrees of freedom T-test. Only 
one grid cell region, with the correlation coefficient of 0.51, is not significantly correlated to the 
solar activity at the 0.05 level of significance; and the remaining ninety-five grid cells are 
significantly correlated to the solar activity at the 0.05 level of significance. 

In summary, the present analysis has demonstrated that T-CPM is significantly sensitive to 
solar activity on the basis of large volume of data. T-CPM changes due to solar activity changes 
taken to bracket the assumed to be associated with the 11 year cycle are large. Our present results 
show that there is a relationship between solar activity and temperature structure in mesopause 
regions. These results demonstrate the potentially significant physical mechanism which may 


link the 11 year solar cycle and the dynamic state of the mesopause. 
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Figure Captions 


Figure 1. Examples of the CPM data quality control process. (a) The temperature profile is 
the 100th event in 75361 orbit at 2015-11-4 night; (b) The temperature profile is the 40th event 
in 76214 orbit at 2015-12-31 night 


Figure 2. Distribution of global available profiles of T-CPM. (a) Amount of total 
temperature profiles of T-CPM for each year before CPM quality control (hollow hexagon) and 
after CPM quality control (solid circle); (b) Distribution of 4,492,664 profiles of T-CPM for each 
integer latitude interval after CPM quality control 


Figure 3. The scatter plot and the annual-mean of T-CPM and solar flux index during 2002- 
2015. (a) The scatter plot of T-CPM and F107; (b) The annual-mean of T-CPM and Fio7 


Figure 4. The annual-mean of H-CPM during 2002-2015 


Figure 5. Latitudinal distribution of solar response and solar cycle dependence with every 
10° interval. (a) Latitudinal distribution of solar response of T-CPM, and the vertical bars show 
the fluctuation range of solar response; (b) Latitudinal distribution of correlation coefficients 


between T-CPM and F107 
Figure 6. Latitudinal distribution of 4,302,764 temperature profiles with every 10° interval 


Figure 7. Latitudinal distribution of 14yr-average for T-CPM and H-CPM with every 10° 
interval. (a) Latitudinal distribution of 14yr average of T-CPM; (b) Latitudinal distribution of 
14yr average of H-CPM 


Figure 8. Maps of the solar response and solar cycle dependence. (a) Spatial distribution of 
the solar response of T-CPM, and contour interval is 1 K/100SFU; (b) Spatial distribution of the 


correlation coefficients between T-CPM and F107, and contour interval is 0.05 above 0.53 


Figure 9. Maps of 14yr-average for T-CPM and H-CPM. (a) Spatial distribution of 14yr 
average of T-CPM, and contour interval is 2K; (b) Spatial distribution of 14yr average of H- 
CPM, and contour interval is 0.5km 
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Figure 1. Examples of the CPM data quality control process. (a) The temperature profile is the 100th event in 75361 


orbit at 2015-11-4 night; (b) The temperature profile is the 40th event in 76214 orbit at 2015-12-31 night 
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Figure 2. Distribution of global available profiles of T-CPM. (a) Amount of total temperature profiles of T-CPM for 
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Figure 3. The scatter plot and the annual-mean of T-CPM and solar flux index during 2002-2015. (a) The scatter plot 
of T-CPM and F07; (b) The annual-mean of T-CPM and Fi 
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Figure 4. The annual-mean of H-CPM during 2002-2015 
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Figure 5. Latitudinal distribution of solar response and solar cycle dependence with every 10° interval. (a) 
Latitudinal distribution of solar response of T-CPM, and the vertical bars show the fluctuation range of solar 
response; (b) Latitudinal distribution of correlation coefficients between T-CPM and Fio7 
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Figure 6. Latitudinal distribution of 4,302,764 temperature profiles with every 10° interval 
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Figure 7. Latitudinal distribution of 14yr-average for T-CPM and H-CPM with every 10° interval. (a) Latitudinal 
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Figure 8. Maps of the solar response and solar cycle dependence. (a) Spatial distribution of the solar response of T- 
CPM, and contour interval is 1 K/100SFU; (b) Spatial distribution of the correlation coefficients between T-CPM 
and F107, and contour interval is 0.05 above 0.53 
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Figure 9. Maps of 14yr-average for T-CPM and H-CPM. (a) Spatial distribution of 14yr average of T-CPM, and 
contour interval is 2K; (b) Spatial distribution of 14yr average of H-CPM, and contour interval is 0.5km 
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